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Abstract 

The boundary supersymmetric sinh-Gordon model is an integrable quantum field 
theory in 1 + 1 dimensions with bulk N = 1 supersymmetry, whose bulk and boundary 
S matrices are not diagonal. We present an exact solution of this model. In particular, 
we derive an exact inversion identity and the corresponding thermodynamic Bethe 
Ansatz equations. We also compute the boundary entropy, and find a rich pattern of 
boundary roaming trajectories corresponding to c < 3/2 superconformal models. 
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1 Introduction 



The supersymmetric sinh-Gordon (SShG) model [0-0 is one of the simplest examples of a 
1 + 1 dimensional integrable quantum field theory with N — 1 supersymmetry. Indeed, the 
particle spectrum consists of one Boson and one Fermion which have equal mass and which 
enjoy factorized scattering J7j. As such, SShG is a valuable toy model. 

In this article, we consider the boundary SShG model, with boundary conditions that 
preserve the bulk integrability, but not necessarily the bulk supersymmetry [ 



111. In 



addition to its usefulness as a simple prototype, we expect that this model may also have 



applications to quantum impurity problems fl2]l . 

An interesting feature of the boundary SShG model is that the S matrices which have 
been conjectured for both bulk and boundary scattering are not diagonal. Our main objective 
is to perform a thermodynamic Bethe Ansatz (TBA) analysis ||, ||13||-||17|| for this model, 
using these conjectured S matrices as inputs. Such analysis can provide checks on the input 
S matrix data, as well as information about the underlying boundary conformal field theory 

Jl- 



Conventional wisdom suggests that the problem of determining the necessary Bethe 
Ansatz equations is intractable, due to the fact that the boundary 5* matrix is not diag- 
onal. Nevertheless, we succeed to determine the Bethe Ansatz equations and carry out the 
TBA analysis for the boundary SShG model. This is the first example of a model with 
non-diagonal boundary S matrix which is exactly solved. 

We also obtain an expression for the boundary entropy p0|,||l7| for the boundary SShG 
model. Moreover, we find a rich pattern of boundary roaming trajectories corresponding 
to c < 3/2 superconformal models [pi ], ||22|| , thereby generalizing previous work on bulk 
23fl-|25l| and boundary p6| , [27] roaming. 

The outline of this article is as follows. In Sec. 2 we review the scattering theory of the 
boundary SShG model, which serves as our input. Here we also show that the strong-weak 
duality symmetry of the bulk model (see, e.g., |28|) extends also to the model with boundary. 
Moreover, we introduce the notations and conventions which are used throughout the paper. 
In Sec. 3 we formulate the so-called Yang matrix [29| and relate it to a commuting transfer 
matrix, which is the true starting point of any TBA analysis. For the problem at hand, 
we require a boundary version of the Yang matrix [[!(J, which presents an interesting 
complication with respect to the more familiar case of periodic boundary conditions. In Sec. 
4 we use the open-chain fusion formula |32| to derive an exact inversion identity, using which 
we obtain the eigenvalues of the transfer matrix in terms of roots of certain Bethe Ansatz 
equations. That such an inversion identity exists is presumably due to the fact that the bulk 
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S matrix satisfies the so-called free Fermion condition 113311 , In Sec. 5 we use these results 



to derive the TBA equations. Certain remarkable identities lead to very simple formulas, in 
particular for the boundary entropy. In Sec. 6 we use our result for the boundary entropy to 
obtain boundary roaming trajectories. Finally, in Sec. 7 we discuss our results and describe 
some possible generalizations. 



2 Review of boundary SShG scattering theory 

In this Section, we review the bulk and boundary S matrices which have been proposed for 
the boundary supersymmetric sinh-Gordon model. As mentioned in the Introduction, these 
S matrices will be used as inputs in the calculations that follow. We also show that the 
strong-weak duality symmetry of the bulk model extends also to the model with boundary. 



2.1 Bulk 



In order to understand the SShG scattering theory, it is essential to first consider a re- 
lated model, namely, the supersymmetric sine-Gordon (SSG) model, whose Euclidean-space 
Lagrangian density is given by 

1 - - M 2 

tfuik = 2 {MM + W + ^) + m ^ cos /^ + ^2 sin2 P<t> > (2- 1 ) 

where 4>(z ,z) is a real scalar field, ip(z , z) and ip(z , z) are the components of a Majorana 
spinor field, and (3 is the dimensionless coupling constant. The Lagrangian density for the su- 
persymmetric sinh-Gordon (SShG) model is obtained by analytic continuation to imaginary 
coupling, i.e. setting (3 — i(3 with (3 real: 

1 - - M 2 

C^shG = (^0gL0 + ipd z tp + ipd,ip) + iMiPi/j cosh (3<p + — - sinh 2 (3<j> . (2.2) 

2 2(3 2 

Both of these models have N — 1 supersymmetry (without topological charge) jl]], [0 and 
are integrable @, [|]. Q 

Observe that SSG has a periodic potential, which admits classical soliton solutions that 
interpolate between neighboring minima. Correspondingly, it has been proposed PS[-P7[ 

^^Of course, a similar relation exists between the usual (non-supersymmetric) sine-Gordon (SG) and sinh- 
Gordon (ShG) models. It is more straightforward to infer the scattering theory for the trigonometric (SG, 
SSG) models than for the hyperbolic (ShG, SShG) models, because the former have kinks with topological 
charge, whose non-diagonal S matrices must satisfy highly restrictive constraints @,|35|,@|. The S matrices 
for the hyperbolic models are inferred by analytic continuation of the corresponding breather S matrices, as 
is explained in more detail below. 
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that the SSG quantum spectrum consists of supersymmetric multiplets of kinks of mass m 
and breathers (bound states of kinks) of mass m n = 2m sm(nan), n = 1 , 2 , . . . , [j-], where 

« = (2-3) 

and [x] denotes integer part of x. Hence, breathers can be present only if < a < |. 
The lightest (n = 1) breathers are identified as the elementary particles (Boson, Fermion) 
corresponding to the fields in the Lagrangian density ( [2.1| ). 

Upon making the analytic continuation to SShG (which is the model of primary interest), 
we see that the potential is no longer periodic, and hence, there are no longer any classical 
soliton solutions. Thus, the SShG quantum spectrum does not contain kinks; it consists only 
of the elementary particles of some mass m corresponding to the fields in the Lagrangian 
density ( |2.2|) , i.e. corresponding to the n = 1 SSG breather. Setting (3 = i(3 in Eq. (|2.3|) , we 
obtain 

t 

a = ^- = -B, (2.4) 

where we have introduced the SShG parameter B. 

Since the SShG spectrum corresponds to the n = 1 SSG breather, we infer that the SShG 
S matrix S(9) for two particles of rapidities 9\ , 9 2 (and corresponding energy = m cosh 0± 
and momentum Pi = msinh^j, i — 1 , 2) is given by the analytic continuation of the n — 1 
SSG breather S matrix, 

S{0) = S ShG (0) S S usy(0) , (2.5) 

where 9 = 9 1 — 9 2 . The scalar factor Sshci^) is given by 

sinh9 -ism(2Bn) 
S S hG{9) = . . . r • (2.6) 

smh 9 + i sm(zi37rj 



This is the S matrix of the usual (non-supersymmetric) sinh-Gordon model [38|,|39], which 



is the analytic continuation of the n = 1 SG breather S matrix ||40|| , but with a different 
dependence on the coupling constant. It satisfies 

S ShG (0) S ShG (-9) = 1 , S ShG (9) = S ShG (™ ~ 0) . (2.7) 

The factor Ssusy{9) is given by f] 

Ssusy(9) = Y(9) R{9) , (2.8) 



2 The matrix Ssusy(9) for SSG was first obtained H in terms of an unknown parameter A by solving the 
constraints coming from supersymmetry and factorization. The identification of A in terms of a was made 
in & 
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where R{9) is a 4 x 4 matrix acting on the tensor product space V <8> V, where V is the 
2-dimensional vector space of 1-particle states. We choose {\b(9)) , \f(9))} to be the basis of 
V (corresponding to a Boson, Fermion with rapidity 9, respectively); and hence, the basis 
of V ® V is given by {|6 a , b 2 ) , \h ,f 2 ),\f l , b 2 ) ,\f\, f 2 )}, where \h , b 2 ) = , b{9 2 )) , etc. 
In this basis, R(9) is given by 



R{9) 



with 



a±{9) = ±1 



2i sin Btt 
sinh 



/ a+{9) d(0) \ 

b c(9) 

c(9) b 

\ d{9) a_(0) / 

% sin i?7r 



(2.9) 



sinh | 



sin i?7r 
cosh | 



(2.10) 



It is important to note that the matrix elements of R{9) satisfy the "free Fermion" condition 

(2.11) 



#3) Rn(9i 



a + a_ + b 2 = c 2 + d 2 . 
The matrix R{9) is a solution of the Yang-Baxter equation Q 

Rn(9i — 9 2) Riz{9\ — 9 3 ) R 2 3(9 2 — 9 3 ) = R 2 2,{9 2 — 9 3 ) R\ 3 {9\ 
This matrix is both P and T invariant, 

V12 Ru(9) V 12 = R 12 {9) , R 12 (9Y^ = R 12 (9) , 
where ti denotes transposition in the i th space, and V is the permutation matrix 



(2.12) 



(2.13) 



V 



( 1 \ 

10 

10 

V 1 / 



(2.14) 



Moreover, Y(9) is a scalar factor given by 



Y(9) 



sinh 



sinh 2 — i sin Bit 



cxp 



dt smh(it9/ir) sinh(t(l + B)) sinh(iB) 



t 



cosh t cosh 2 I 



, (2.15) 



We use the very useful convention (which is standard in the spin- chain literature @ - Q, but 
unfortunately not in the field-theory literature), whereby Rij(0) acts nontrivially on the i th and j th 
vector spaces. For instance, in the Yang-Baxter equation, the R matrices act on V® 3 , and therefore 
Ri2(0) = R{9) ® I , #23(0) = I O R{0), etc., where I is the unit matrix. 
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which is a solution of the unitarity and crossing constraints 



sinh 2 



Y(6) Y{-6) = ^ Y{6) = Y{vk-&). (2.16) 
sinh 2 + sm -o 71 " 

Let us denote the total scalar factor by Z(9) 

Z{9) = S ShG (9) Y{6) . (2.17) 

One can show that Z{6) has the integral representation [H] 

Z ( e) = exp ( r - ^/^inl^l-^sin^A , (2.18) 

sinh 2 + i sin Bn \J t cosh t cosh ^ J 

which is the same as the expression in Eq. (|2.15|) , except with B — > —B. It has no poles [] 
in the physical strip (0 < Im 9 < tt), provided B lies in the range 

0<B<1, (2.19) 

which corresponds to < (3 2 < oo. 

In short, the proposed SShG bulk S matrix S(9) is given by 

S(6) = Z{6) R{6) , (2.20) 



where the scalar factor Z(8) is given by Eq. (|2.18 ) and the matrix R(6) is given by Eqs. 

(PD, O- 



It is known (see, e.g., ^B|) that the SShG bulk S matrix is invariant under the strong- 
weak duality transformation (3 — > 47r//3, which implies 

B^l-B. (2.21) 

Indeed, this invariance can be checked by inspection of the matrix elements ( |2.1(J| ) of R{6) and 
the expression ( p,18| ) for Z{9). (The factors Sshc{9) an d are not separately invariant.) 



Note that this transformation maps the range (|2.19| ) into itself. 
2.2 Boundary 

We turn now to boundary conditions and boundary scattering, following the framework 
developed by Ghoshal and Zamolodchikov H. An investigation of which boundary terms 



'Although Y(9) has a pole at 9 = i2Bir, it is canceled by a corresponding zero of Sshg(6)- 
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can be added to the bulk SShG model (|2.2|) without spoiling (classical) integrability has led 
to the following results ||: the boundary Lagrangian 



LlTndary = A CO S h/3(0 - O ) +M# + ^ + ^- ) M ^ ±1 , (2.22) 

breaks supersymmetry but preserves integrability; and 

JwU = T§ cosh/30 ± ^ (2.23) 







preserves both supersymmetry and integrability. Notice that the boundary terms fl2.22| ) 



involve a total of 5 boundary parameters A , , M , e , e. If e , e are nonzero, then Fermion 
number is not conserved. 

The proposed boundary S matrix S(9) for a particle of rapidity 9 is given by (compare 
with Eq. (p|)) 

s(0) = s ShG (e ; n , •&) sf USY (e ; ^) . (2.24) 

The scalar factor SshG(9 ] V , which depends on two boundary parameters 77,1?, is given 
by 

Ssk G (e ; V ,0) = x (9) x^e ■ -L-) x t (9 ■ ) , (2.25) 

7T 7T 



where 



X (9) = (1)(1 + 2B)(2 - 2B) , X 1 (9 ; F) = _ 1 - , (2.26) 



with 



sinh(, + 

1 ; sinh (| - if) K } 

This is the boundary S matrix of the usual (non-supersymmetric) boundary sinh-Gordon 
model, which is the analytic continuation of the n — 1 boundary sine-Gordon breather S 



matrix 44 . It satisfies 



?7T ?7T 

S ShG (9 -,71,0) S ShG (-9 ; V , 0) = 1 , S s/lG (- + 0; , 0) S 5/lG (2#) = S SftG (y -0;t],&). 



(2.28) 



5 We make an effort to distinguish boundary quantities from the corresponding bulk quantities by using 
sans serif letters to denote the former, and Roman letters to denote the latter. 

6 The boundary S matrix is obtained in terms of a set of boundary parameters [rj , 1? , ip , s) by solving the 



boundary Yang-Baxter equation. The relation of these parameters to those in Lf^^ (2.22) is not known 
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The factor S^ SY (9 ; (f) is given by |10 



Sf USY {9^)=^\9- v ) R^(9;<f), (2.29) 

where e is a discrete parameter which can be either +1 or — 1, and if is a continuous boundary 
parameter. R^ e ^ (9;(f) is a 2 x 2 matrix acting on the vector space V of 1-particle states, 
which is given by 

I cosh § + i sinh | -si sinh , o 



where 



-i sinh cosh f Gf - z sinh | G[ 



r f sinh o? + ee £lf --. smh . 2 D ^ if e' = —e 

V ' l+esiniJ7r J 



and 



2(e + sinfivr) 



(2.32) 



sin Bix 

The matrix R^ £ \9 ; </?) is a solution of the boundary Yang-Baxter equation j|5 

#i2(0i - 2 ) RS £ Vi ; ¥>) #12(^1 + e 2 ) R<i\o 2 ; a.) 

= R ( 2 £ \e 2 ; 09) i2 12 (0! + e 2 ) r{ £] (9 1 ■ <p) R^e, - e 2 ) . (2.33) 

We remark that for 09 — > ±00, the matrix R( £ \6;(p) becomes diagonal and commutes with 
linear combinations Q ± Q of the supersymmetry charges |TH|, ||11| . 

In order to determine the scalar factor Y^ £ \9 509), we recall that the full boundary S 
matrix must satisfy boundary unitarity S(9) S(—9) = I and boundary cross- unitarity |§, 
which can be written in matrix form as 

tr S (y + 6)^ V01 S 01 (29)^ = Si(y - 9) . (2.34) 

We observe that the matrix R^ £ \9 ; 09) satisfies 

R {£) (9 ; <f) R {£ \-9 ; V ) = h{9) I , (2.35) 

where 

h{9) = ico + d sinh 2 - + c 2 sinh 4 - J cosh 9 , (2.36) 
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and 



, r 2 COsh 2 ip if S = +1 r 2 e £2 V r 2 g £2^ 

C °H 2 • u2 . c n > C! = — + 2e, c 2 = — . . 2.37 

r 2 smh </? ife = -l l + £SinS7r (1 + e sin I?7r) 2 



Also, 



where 



Setting 



tr R? } (y + 9 ; </?)*° P i ^oi^)' 1 = 9(0) ^(y -*;¥>), (2-38) 



£ sinh — « sin 7r£> , 

9 < 9 > = ' (2 ' 39) 



Y<<>(9 ;<,)=¥<"(«) Y<"(e ;¥ .), (2.40) 
it follows that (9) and (9 ; </?) must satisfy 

Yf Y<?>(-0) cosh^ = 1 , Y<?>(y + 0) F(20) </(*) = Y^ - 9) , (2.41) 

and 

'( e )/7).,„Yv( e )f_/3 ■ ,^ _l „. 0^1,2 £ , „_ ™u4 <M _ i v( £ )/'!![ _l a.,^ - v( £ )^ 77 



(0 ; Yj" (-9 ;<p) ( c + Cl sinh 2 - + c 2 sinh 4 -1=1, Yj" (y + fl;^=Yl e ^-J; V ) ) 



(2.42) 



respectively. 

For simplicity, we shall henceforth restrict our attention to the case e — +1, and so we 
shall drop the superscript (e). We propose the following integral representations for Y o (0) 
and Y x (#;</?): 



Yo(0) = ^ . , /fl . -T exp 



1 Z" 00 sinh(2^/7r) sinh(t(l + B)) smh(tB) 



v^sinhd + f) \2J t cosh 2 tcosh 2 | 

= 1 ( r°° dt cosh(tC/vr) cosh(|(l - 2B)) sinh(|(l + f )) sinh(f ) 

r cosh \ Jo t sinh t cosh | 

(2.43) 

where C is a function of the boundary parameter p> defined by 

( = cos' 1 (1 + e" 2v (l + sin^vr)) . (2.44) 
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In order to streamline the notation, let us denote the set of boundary parameters {i] ,ip} 
by £, and denote the total scalar factor by Z(6 ; £) 



Z(0;O = S Sha (e;ri,#) Y(6;<p) 

= Z (6)Z 1 (9;{). (2.45) 

The proposed SShG boundary S matrix is then given by 

S(0;£) = Z(0;O R(9;<p). (2.46) 

We now observe that the boundary S matrix is also invariant under the strong-weak 
duality transformation ( 2.21|) . Indeed, it is evident that the matrix R(8 ; if) ( |2.30| ) has this 



invariance, if we assume that the parameter (p remains invariant under this transformation. 
Let us now consider the scalar factor. The part of the scalar factor that does not depend on 
boundary parameters can be written in the form 

Z o (0) = MO) MO) (2-47) 

i 

V2sinh(| - f ) 

x exp f _i r ± ^H^o/n) [cosh(t(1 _ 2£?))(1 + 2cosht) + coshi] w 48) 

\ 4 J t cosh t cosh | I 

in which the duality invariance is manifest. (The factors Xo(^) and Yo(6 l ) are not separately 
invariant.) Finally, the part of the scalar factor which does depend on boundary parameters, 

AnB AifiB 

z 1 (6 -,0 = MO; -H MO ; ) Mo ; <p) , (2.49) 

7T 7T 

is also invariant under duality, since each of its factors are separately invariant (provided the 
boundary parameters rjB and dB are assumed to remain invariant). 



3 Yang equations 

Having specified the bulk and boundary S matrices, we are ready to start the TBA program. 
The first step is to formulate the Yang matrix and relate it to a commuting transfer matrix. 
Since this is not obvious for the case of boundaries, we begin by reviewing the more familiar 
case of periodic boundary conditions. 
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3.1 Closed 



We consider N particles of mass m with real rapidities Q\ , . . . , On and two-particle S matrix 
S(0) (|2.20|) in a periodic box of length L >> — . The Yang equation |29|1 , Hl6| for particle 1 
(which has momentum Pi = msinh#i) is given by 



( 



ALm sinh ( 



^Y {l) -i)\e 1 ,...e N ) = Q. 



where Ym is the "Yang matrix"[] 

= Sin(0i — On) <Si,Ar-i(0i — N -i) • • • £12 (#i — #2) 



(3.1) 



(3.2) 



which acts on V® N . There are similar equations, and corresponding matrices Yu\, for the 
other particles % — 2 , 3 , . . . ,N. 

The objective is to diagonalize Yuy The key to this problem is to relate Yu\ to an inho- 
mogeneous closed-chain transfer matrix, for which there are well-developed diagonalization 
techniques. (For reviews, see e.g. fy] - f43|.) Indeed, consider the transfer matrix (see 
Figure |I|) 





















) 




• • • 







N 



2 1 

Figure 1: Closed-chain transfer matrix. 



Tdosed{®\Q\ , On) = tr {5'oiv(6 l — 8 N ) ■ ■ ■ S 02 (9 — 6 2 ) S i(6 — Ox)} , (3.3) 

with inhomogeneities 0% , ■ ■ ■ , On- Notice that we have introduced an additional ("auxiliary") 
2-dimensional vector space denoted by 0. The product of 5* matrices inside the trace (the so- 
called monodromy matrix) acts on \/®( Ar + 1 ); but after performing the trace over the auxiliary 
space, one is left with an operator which acts on the ("quantum") space V® N . Because S(0) 
satisfies the Yang-Baxter equation, the transfer matrix commutes for different values of 

[^closed 

(0\0i , • • • , N ) > ^closed 

{0'\0 X ,... ,0 N )] =0. (3.4) 

7 We remind the reader that we are using the convention explained in Footnote [| 
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Let us now evaluate this transfer matrix at 9 = Q\. Using the fact that 5(0) = V (the 
permutation matrix ( p.14 )) and V 2 = I, we see that 



Tdosed^ ,...,9 N )= tro{(Poi^Ol)5oAr(^ - 9 N ) ■ ■ ■ (V lV l) S 02 (9 1 - 9 2 ) V l} • (3.5) 

Finally, using V 01 S 0i V 01 = S u and tr P i = h, we conclude that 7- dosed (0i|6>i , . . . ,9 N ) = 
Km. In general, we have 

K (i) = T doaed (0 i \O 1 , . . . ,9 N ), i = l,...,N. (3.6) 

This is the sought-after relation. In order to diagonalize the Yang matrices Y(a, it suf- 
fices to diagonalize the commuting closed-chain transfer matrix T c i ose d(9\9i ,.. . ,9 N ). That 
calculation, as well as the corresponding bulk TBA analysis, is described in ||. 

3.2 Open 

We now turn to the case with boundaries, which is our primary interest in this paper. We 
therefore consider iV particles of mass m with real rapidities 9i , . . . , 9^ in an interval of 
length L » i with bulk S matrix S{9) (lZ20|) and boundary S matrix S(0 ; (|^46|) . The 
Yang equation for particle 1 is given by |3(| , pT| 

[e 2iLmsM F (1) — I) \9i , . ..9 N ) — 0, (3.7) 

where the Yang matrix Km is now given by 

K (1) = S^x ; e_) 5 2 i(0! + 6 2 ) ■ ■ ■ S N1 (9 1 + 9 N ) S 1 (9 1 ; ^+) SWfa " M ' • ■ S 12 {6 1 - 9 2 ) ,(3.8) 

where the subscripts ± denote the left and right boundaries. (There are similar matrices 
Kj) for the other particles.) In analogy with the case of periodic boundary conditions, the 
key to diagonalizing the Yang matrix is to relate it to an inhomogeneous open-chain transfer 
matrix |46| (see Figure 0) 




N 2 1 

Figure 2: Open-chain transfer matrix. 
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T{e\e 1: ... ,e N ) = tr {s M + m;£ + )' s 0N (e - e N ) ■ ■ ■ s Q1 (e - e 1 ) 

x s o (0;£_) s 01 (e + e 1 )---s 0N (e + e N )}, (3.9) 

which commutes for different values of 9 

[t(6\0 u ...,0 n ),t(&\6 1 ,...,6 n )] = Q. (3.10) 



Using the boundary cross-unitarity relation (|2.34| ) as well as the Yang-Baxter equation ( 2.12|) , 
( |2.33| ), one can show that 

Y (f) =T(0 i \0 1 ,... ,9 N ), i = l,...,N. (3.11) 

A proof for the case N = 2 is presented in Appendix A. Q Hence, in order to diagonalize 
the Yang matrices Yu\, it suffices to diagonalize the commuting open-chain transfer matrix 
t(9\9i , . . . , 9n). It is to this task that we devote the following section. 



4 Inversion identity and transfer- matrix eigenvalues 

In this Section, we consider the problem of determining the eigenvalues of the inhomogeneous 
open-chain transfer matrix ( |3.9|) . Our approach will be to first derive an exact so-called 
inversion identity. This approach has been used in the past to diagonalize simple (e.g., 
Ising) closed-chain transfer matrices [^7|,[15|,@. 



4.1 Inversion identity 

Instead of working with the "dressed" transfer matrix ( |3.9|) , it is more convenient (see Foot- 
note |^ below) to strip away the scalar factors from the bulk and boundary S matrices, and 
to work instead with the "bare" transfer matrix 

t(0|0i,... ,9 N ) = tro{R (-# + 27r;^ + ) to R m {9 - 9 N ) ■ ■ ■ R m {9 - 9{) 

x R o (0 ;</?-) Roi(0 + 9 1 )---R ON (9 + 9 N )}, (4.1) 

where R{9) is given by ( |2.9|) and R(9 ; ip) is given by (|2.30| ) with e = +1. 



There are two key points involved in obtaining the inversion identity. The first key point 
is to observe that the bulk S matrix degenerates into a one-dimensional projector for a 



8 We therefore fill a gap left open in p0|| , where it was first observed that the open-chain Yang matrix is 
related to the Sklyanin transfer matrix; but neither the precise form of the relation nor its proof was given. 
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certain value of 9 (= — iir): 



S(—iir) oc 





/ 


1 








-1 \ 


1 
















2 


















V 


-1 








1 / 



(4.2) 



Hence, it is possible to "fuse" ||48fl, |[49|| ,||41|| in the auxiliary space, and thereby obtain a fusion 
formula of the form 



t(0|0i 



, e N ) t(e + m\e x , . . . , e N ) a i{e\9 x , . . . , e N ) + a 



(4.3) 



where \{6\6i , . . . ,9^) is a "fused" open-chain transfer matrix (see Figure |3|), and here A 
represents a product of certain quantum determinants [50fl. The fused transfer matrix is 
constructed from the fused bulk S matrix R{9) and the fused boundary S matrix R{9 ; ip) , 
using the "fused" 3-dimensional (instead of 2-dimensional) auxiliary space. 




Figure 3: Fused open-chain transfer matrix. 

The second key point is that both R(6) and R(9 ; cp) can be brought to upper-triangular 
form by a ^-independent similarity transformation. This remarkable fact is presumably due 
to the fact that R(9) satisfies the free Fermion condition ( |2.11| ) (cf, |f£3}|). As a result, the 
fused transfer matrix is proportional to the identity matrix 

t(0|0i,... ,9 N ) ocl. (4.4) 

It follows from the fusion formula that the transfer matrix obeys an exact inversion identity 

t(0|0! , . . . , 9 N ) 1(9 + Z7r|0 a , . . . , 9 N ) = f(9)l , (4.5) 

where f(9) is a calculable scalar function. We find (see Appendix B for more details) 



/(*) 



16 sinh 2 9 



sinh(0 - iBn) sinh(0 + iBn){l + sin Bit) 2 sin 2 Bit 
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.7 = 1 



^ cosh(§(0 - 6j) - iBn) cosh(±(0 + 9,) - iBn) 
cosh(±(0 - Bj)) cosh(±(0 + 0j)) 



-j^j- sinh(i(0 - 0,-) - z^tt) sinh(±(0 + 9 3 ) - iBir)\ 
' 11 rinh(J(e-^)) ^(§(0 + 0,)) J 

Wfl 1 9 j )+iB 7 r) 



7 = 1 



- 0j) + zBtT 

1 cosh(±(0 + 0j 



) COsll(; 



/ Jt cosh(|(g- 

1 1 cosh(|(^-%)) 

sinh(i 



+ vLBtt) 



where 



a 



I ? 7T 

= sinh 2 (-(M + y -0)) 
= sinh 2 (-H7r5 + y + tf)) 



c = sinh 2 (-(m5 + — + 6)) 
2 2 

= smh\hinB-^ + 9)) 



sinhd^ - //' ' 

I <9<7r I o -yr 

sinh 2 (-(m J B + -)) + e^" S inh 2 (-(-wB + y + 0)) x 

! ZTT I Z7T 

smh. 2 (-(inB + — )) + sinh^-^TrS + 9)) x [<p 

2 2 2 2 

<^sinh 2 (i(M-y + 0)) x [(p_ 



-<p- 



e- v - sinh 2 (-(nr J B + 

V 2 V 2 



V 2 V 



1/. ^ «7T N 



lit . , 

sinh' 



^ /iwl r 

y + ^)) x [p_ 



Notice that the function f(9) is invariant under the duality transformation B — > 1 — 5. This 
inversion identity is one of the main results of this paper. We have checked it numerically 
up to N = 3. 



4.2 Eigenvalues 

We now proceed to determine the eigenvalues of the transfer matrix. First, observe that 
by virtue of the commutativity property ( |3.10|) , the bare transfer matrix t(0|0j , . . . , 0jv) has 
eigenstates \9\ , . . .9n) which are independent of 9, 

1(919, ,9 N )\9 1 ,...9 N ) = £(0)0! , . . . , 9 N )\9, , ... 9 N ) , (4.8) 

where £(0|0i , . . . , 0jv) are the corresponding eigenvalues. Acting on |0i , . . . 9^) with the 
inversion identity, we obtain the corresponding identity for the eigenvalues 

£(0|0x , . . . , 9 N ) 2(9 + m\9, ,...,9 N ) = f(9) . (4.9) 



14 



Moreover, one can show that the bare transfer matrix t(6\6i , . . . , 9jsr) is a periodic func- 
tion of 9 with period 2ni f\ 



t(0 + 2iri\6i,... ,9 



N 



Wi 



whose asymptotic behavior for large 9 is given by 



where 



,9 N )~-e 3e I for 9^oo 



(1 + sin Bti) sin Bn 
Correspondingly, the eigenvalues obey 



(4.10) 



(4.11) 



(4.12) 



£(0 + 2tti|0i , . . . ,9 N ) 
2(9\9 U ... ,9 N ) 



■ 9n) , 



~ — e se for 9^oo 
32 



(4.13) 



The eigenvalues £(0|0i , . . . , 9 N ) are uniquely determined by the zeros and poles of f(9), 
together with periodicity and asymptotic behavior. Indeed, observe that f{9) is a product 
of two factors. Let 9 = z£ , z^ be zeros of the first, second factors, respectively. Then z£ 
obeys 



N 

n 



tanh(±(4 - 6j) - iBit) tanh(±(z+ + 9j) - iBti) sinh 2 (±(m£ + f - «+)) 



tanh(|(4 - S )) 



tanh(|(4 + %)) 



sinh 



2/1, 



-M + f + 4)) 



x 



e~f- sinh 2 (|(M + f )) + e^- smh 2 (|(-M + f + z+)) 
e-f- smli 2 ^^ + f )) + e*- sinh 2 (±(«7r.B + f - 4)) 



x [v9_ - ^ + ] ,(4.14) 



and z fe obeys 



n 

i=i 



tanh(|(Zfc - 0j) + iBvr) tanh(±(^ + 9 3 ) + iBvr) sinh 2 (i(m5 + f + z")) 



tanh(±( % 7-0,)) 



tanh(|(^ + %)) 



sinh 2 (i(z7rfl-f + 



e - * 1 - sinh^KivrS + f )) + e<^ sinlT(§(i7r.B - f + z;:)) 
e~^- sinh 2 (|(z7rfi + f )) + ~&- sinh 2 (±(M + f + ^j) 

These are our "magnonic" Bethe Ansatz equations. 



x [y>_ - <p + ] . (4.15) 



9 This is not the case for the dressed transfer matrix r 
factors. 



,9 n), due to the presence of the scalar 
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It follows [J] that f\6) can be represented as 

f(Q) = -— sinh 2 9 ^=° Sinh ^ - Z ^ Sinh ^ + Z ^ Sinh ^ - ^ Sinh ^ + ^ (4 16) 
' 16 Yl^ =1 smh 2 (6 -6 k )smh 2 (6 + 6 k ) 

It now follows by similar arguments that 

me e ) = csinhg H^^iMg - sinh ^ 9 ± z ^ sinh ^ g - » sinh (^ + <jH 

nti sinh(i(# - 4)) sinh(±(0 + 6 k )) cosh(±(fl - cosh(±(# + 6 k )) ' 

(4 

is the unique solution to the inversion identity (|4.9| ) with the properties (|4.13|) . Note that 
there are N + 1 pairs of roots zjf , whereas in the case of periodic boundary conditions || 
there are only N. The appearance of the additional pair of roots is due to the fact that 
the boundary S matrix R(6 ; ip) is not diagonal. The existence of these roots is essential for 
obtaining the correct asymptotic behavior; and it can be easily checked for the case N = 0. 

In summary, the eigenvalues of the bare transfer matrix (|4.1| ) are given by Q4.171 ), where 
4 satisfy Eqs. fl£I3D , (|4TT5| ). 



4.3 Structure of Bethe Ansatz roots 



Before performing the thermodynamic (N — > oo) limit (which is the subject of the next 
section), it is necessary to first understand the structure of the Bethe Ansatz roots. Following 
0, we observe that the Bethe Ansatz Eqs. ( f4.14|) , (|4.15 ) have roots of the form 



x k + iBir 
x k + iBn + m 



z, 



x k — iBir 
x k — iBn — m 



(4.18) 



where x k are real and satisfy 

x r tanh(|(xjt — 8j — iBn)) tanh(|(xfe + 6 



n 

i=i 



j 



iBn)) 



tanh(|(a;fc — 9j + iBti)) tanh(|(xfc + 9j + iBir)) 



sinh 2 (|(f 



+ x k )) 



sinh 2 (i(f 



e-f- sinh 2 (±(^5 



)) + e ip ~ sinh 



2( I fin 



Xk)) 



e-f- sinh 2 (f (M + f )) + e v- sinh 2 (|(f + x k )) 



-x h )) 
x [<P- -> ¥>+] = 



A- 



0,1, 



AT. 



(4.19) 

10 Indeed, let us denote the right-hand-side of Eq. (|4.16|) by We observe that both /(0) and F(0) 

have the same periodicity (namely, in, which is half the period of £(0|0i , ■ ■ ■ , On)), the same zeros and poles 



in the strip 



< < 



and the same asymptotic behavior. (The apparent poles of /(0) at = ±H3tt 



are canceled by corresponding zeros.) Hence, the function g{9) = F(9)/f(9) is regular everywhere in the 
complex 9 plane, and thus must be constant by Liouville's theorem. By considering the limit 9 — > oo, we see 
that this constant must be 1. 
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Evidently, for each x k , there are 4 possible combinations of roots (z k ,z k ). However, by 
considering the limit B —>■ 0, one can argue that only 2 of these combinations are allowed, 
which we denote by = +1 and e k = —1, respectively: 



e k = +1 : (zl = x k + iBn , z k — x k — iBn - in) , 
e k = -1: (zt = x k + iBn + in , zl = x k - iBn) . 



(4.20) 



Hence, the eigenvalues are specified by {x k ,e k }, k = , . . . ,N: 



Yl k= iI si M0-0 k )smh(9 + 9 k 



z{9\o x ,... ,e N ) 



eo--e N 



csinh#- 



where 



\ e (6) = sinh(^(# - tiBn)) cosh(~(0 + eiBn)) 



(4.21) 



(4.22) 



e k — ±1, and x k satisfy ( |4.19| ). 

To close this section, we observe that the "dressed" transfer matrix ( |3.9|) is simply related 
to the "bare" transfer matrix (|4.1|) by 



JV 



TWW X 



, 9 N ) = Z(-6 + in ■ £ + )Z(6 ; £_) J] [Z(6 - k )Z(6 + 6 k )} l(0|0 a ,...,6 N ), (4.23) 



k=l 



where the scalar factors Z(6) and Z(6*;£) are introduced in Eqs. ( |2.17| ), fl2.45| ). Hence, the 
eigenvalues A(^|^ x , . . . , 6^) of r(9\9i , . . . , 6^) are given by 



N 



Am,...,M £ra = A °^ +) Z ^-l Z{ - 2e) \{z{e-e k )z{e + e k ) 



x £(8\8i , . . . , 6 N ) €Q ... eN , 

where £(9\9i , . . . , 6 l Ar) eo ... ejv is given by (|4.21| ). Here we have used the fact 

Z(6;OZ(-29) 



9(f~ 



(4.24) 



(4.25) 



which follows from the bulk unitarity ( [2.7|) , ( |2.16|) and boundary cross-unitarity (|2.28|) , (|2.41| ), 
( p^ ) relations. 



5 Thermodynamic Bethe Ansatz analysis 

Having obtained the eigenvalues of the transfer matrix and the Bethe Ansatz equations, we 
can proceed to the derivation of the TBA equations and boundary entropy. We begin by 



17 



briefly reviewing the general framework. Following ||,[r7 we consider the partition function 
Z + _ of the system on a cylinder of length L and circumference R with left / right boundary 
conditions denoted by ± (see Figure f|) 



R 




Figure 4: Cylinder on which the partition function Z + _ is defined. 

Z + _ = tre-^-^e-^ 
= {B + \e- LHp \B_} 

w (B + \0)(0\B_}e- LEo forL^oo. (5.1) 

In the first line, Euclidean time evolves along the circumference of the cylinder, and H + ^ is 
the Hamiltonian for the system with spatial boundary conditions ±. In passing to the second 
line, we rotate the picture, so that time evolves parallel to the axis of the cylinder; Hp is 
the Hamiltonian for the system with periodic boundary conditions, and \B±) are boundary 
states which encode initial/final (temporal) conditions. In the third line, we consider the 
limit L — > oo; the state |0) is the ground state of Hp, and Eq is the corresponding eigenvalue. 
The quantity ln(.B + |0)(0|£>_) is the sought-after boundary entropy pO|],[pJ|. Taking the 
logarithm of the above expressions for the partition function, one obtains 

- RF^ -LE Q + ]xi(B + \0) (0\B-). (5.2) 

Whereas the free energy F has a leading contribution which is of order L, here we seek the 
subleading correction which is of order 1. 



n More precisely, we shall compute the dependence of the boundary entropy on the boundary parameters. 
The term in the boundary entropy which is "constant" (independent of boundary parameters) seems to be 
difficult to compute even for simpler models ff7|| , [j5f| . 
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5.1 Thermodynamic limit 



We proceed to compute F using the TBA approach 0, [|Kj-[|17|]. To this end, we introduce 
the densities P±{9) of "magnons", i.e., of real Bethe Ansatz roots {xk} with = ±1, 
respectively; and also the densities pi{9) and p{9) of particles {9 k} and holes, respectively. 
Computing the logarithmic derivative of the "magnonic" Bethe Ansatz equations ( 4.19Q , we 
obtain 12 



P + (9) + P49) 



+ 



— / d& px(&) [$(9 - 9') + $(0 + 9')] 
^ Jo 

-L [-m + 29{6) + ^ v+ (9) + %_(9)] 



(5.3) 



where 



1 d , /tanh(i(0 - iS7r)) 
In 



In — — = 

idO \t&nh(±(9 + iBir)) J cosh 20 - cos 2Bn 



4 cosh 9 sin i?7r 



1 9 



smh(|(| + 0)) 
sinh(i(f -(9)) 



1 



cosh0 



ld_ 
~id9 



In 



e-* 5 sinh 2 (|(z7rS + f )) + sinh 2 (±(f - 0)) 
e-^ sinh 2 (i(«7r5 + f )) + sinh 2 (±(f + 0)) 



4 cosh 6* cos C 
cosh 29 + cos 2C 



(5.4) 



where ( is defined in ( J2.44 ) . Defining pi{9) for negative values of 9 to be equal to pi(|0|), we 
obtain the final form 



P + (9) + P_(0) = i ( Pl * $) ((9) + -L [-$(0) + 2¥(0) + v^ + (0) + VW] 



where * denotes convolution 



d9' f{9 - 9')g{9') . 



(5.5) 



(5.6) 



We next consider the Yang equations, which imply (see Eqs. 



,(ED) 



J2iLm sinh 61 



A(0fc|#l , . . . , 0jv) — 1 



fc = 1 iV. 



where A(0|0i , . . . , 9n) is the eigenvalue of the dressed transfer matrix t{9\9\ , . 
is given by ( |4.24| ). Computing the logarithmic derivative, we obtain 



(5.7) 



9n), which 



Pl (9)+p(9) = ^-\2mcosh9+ [ d9' p x {9') [§ z {9 - 9') + § z {9 + 9' 



)] 



12 The term — 2^l${9) originates from the exclusion of the Bethe Ansatz root Xk = 0. 
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+ 



P+(9')<S> + (9 - 9') + P_(0')$_(0 - 6') 



+ P_(0')$+(0 + 9') + P+{9')$-{9 + 6') 
1 



+ 



-® z {9) - 2$ z (29) + —Im In Z(0 ; £+) + -^Im In Z(0 ; £_) 



(5.1 



where 



§ z {9) = —ImlnZ(9), $±(0) = — ImlnA±(0) 



(5.9) 



Using the fact $±(0) = ±|$(0), and defining P±{6) for negative values of to be equal to 
~~ -P±(|0|)) we obtain 



777 ! 1 

Pi(0)+/3(0) = -cosh0 + — ( Pl * $ z ) (0) + — ((P + -P_)*d>)(0) 

7T Z7T 47T 



+ 



1 



2ttL 



-<$> z {9) - 2<$> z (29) + — Im In Z(0 ; £ + ) + ^™ ln ! ^ 



;5..io) 



We now use Q5.5|) to eliminate P_, and use the expressions (|2.45|) , Q2.49|) to separate the 
various factors in Z(0 ; £) to obtain 



Pi(O)+p(0) 



m 



71 



1 



1 



cosh 9 H P + * $ H pi * $ z 



2tt 



2tt' 



4?t 



+ 



+ 



+ 



1 



2ttL 
9 



— I $2 $ * $ 

4tt 



^/mlnZ o (0) - $> z {29) - ^ * $ 



50 

1 d_ 
~id9 



Im\nY l (9; ( p+) -— *$ 



5 



JmlnYi(0; — * $ 

4tt 



71 



71 



i 89 



71 



7T 



(5.11) 



Noting the "bulk" identity |27 



$z(8) ($*$) (0) = 

47T 



and its boundary counterparts 



A Jm l nYl (0;^-i-(^*$)(0) = 0, 
^/mlnZ o (0)-$z(20)-^(vI/*<|.)(0) = -i 



(5.12) 



$(0) + tf(0) 



(5.13) 
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we remain with the rather simple result 

777 ] j 

pM + W) = -coshe + — (p + *®)(6) + -— 

7r lix ZnL 



+ k{ 



+ k(0\ , 



4r?_P, 



+ K 



7T 



7T 



7T 



7T 



where 



K{9 ' F) ~ z^ lnXl( ^' F) " cosh 2£ + cos ttP 



(5.14) 



(5.15) 



The thermodynamic limit of the magnonic Bethe Ansatz equations and the Yang equa- 
tions, given by (|5.5|) and ( |5.14|) , respectively, are the main results of this subsection. Notice 
that the former depends on the boundary parameters (p±, while the latter depends on the 
(boundary sinh-Gordon) boundary parameters r]±,$±. 



5.2 TBA equations and boundary entropy 

The free energy F is given by 

F = E-TS, 
where the temperature is T = ~, the energy E is 

N 



E = ^^mcosh^fc = — / d9 pi(9)m cosh. 9 . 
k=i 2 



and the entropy S is |13|,|jT5| 

S = —J d9 + p)ln(pi +p) -pilnpi -plnp 

+ (P+ + P_) ln(P + + P_) - P + In P + - P_ In P_| . 
Extremizing the free energy (8F = 0) subject to the constraints 



5P_ = SP + + —5 Pl * $ , 

27T 

5p = -5pi + —SP + * $ , 

27T 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



(which follow from Eqs. ( |5.5| ), ( |5.14| ), respectively) we obtain a set of TBA equations which 

)unda 



is the same as for the case of periodic boundary conditions [pl,pv| 

1 

1 



r cosh 9 




e 2 (9) 



2n 



($ * L 2 ) (9) , 
($ * Li) (0) , 



(5.20) 
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where 



U{9) 



In (l + e~ e ' (e) ) , r = mR. 



Pi 



£2 



in I | 



(5.21) 



We next evaluate F using also the constraints ( |5.5| ), ( |5.14| ) and the TBA equations. From 
the boundary (order 1) contribution, we obtain (see Eq. ( |5.2j )) the boundary entropy 



ln(fl + |0)(0|B_) = — 



de 



+ 2¥(0) + k(0 ; ^) + k( 

2 7T 



71 



+ K 



+ k( 



71 



71 



HQ 



+ [-$(0) + 2tf(0) + ^ + (0) + tf v _(0)]L 2 (0). (5.22) 

In particular, the dependence of the boundary entropy of a single boundary on the boundary 
parameters is given by 



1 

Arc 



dO 



k{9;^-)+k( 



4i$B. 



71 



71 



I*(9) + %(0) L 2 {9) 



(5.23) 



where the kernels k(9;F) and ^(0) are defined in Eqs. (|5.15|) and (|5.4j) , respectively. The 
term involving Li, which had previously been conjectured [p7fl , depends on the boundary 
sinh-Gordon parameters 77,1?. The term involving L2, which had not been anticipated, 
depends on the boundary parameter ip (which appears in R(9 ; ip), i.e., the non-diagonal part 
of the boundary S matrix). This expression for the boundary entropy is another of the main 
results of this paper. 



6 Boundary roaming trajectories 

One application of our result ( |5.23| ) for the boundary entropy is to obtain boundary roaming 
trajectories corresponding to c < 3/2 superconformal models. In order to best explain this 
result, it is helpful to first recall earlier work on bulk and boundary roaming. 

Al. Zamolodchikov [23j first considered the TBA equations for the bulk ShG (non- 
supersymmetric) model with the coupling constant 7 analytically continued to complex val- 
ues, 

7=|±i0o, 0o »1- (6.1) 



13 For the case e = — 1, we obtain a similar result, except the parameter £ appearing in the kernel ^ v (0) 



is now given by C = cos" 1 (-1 + e 2 f(l - sin Bit)) instead of by Eq. fl2.44|) 
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The corresponding effective central charge c e ff(r) interpolates ("roams") between the values 

6 



1 



V 



3,4,5 



(6.2) 



p(p + l) 

corresponding to the unitary c < 1 minimal models [[18]. Indeed, a plot of c e //(r) vs. log(r/2) 
reveals a "staircase" with plateaus at values of c e //(r) equal to c p . 

This result was later generalized |26] to the boundary ShG model: choosing the value of 
r so that c e //(r) lies on some plateau, the boundary entropy sb(F) (where F is a bound- 
ary parameter) interpolates between values corresponding to various conformal boundary 
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The original work |23| was also generalized |27| to the bulk SShG (supersymmetric) 
model. The TBA equations with a similar analytic continuation of the coupling constant 



7T 



B 



71 



±l9r 



» 1 



cause the effective central charge c e //(r) to interpolate between the values 



1 - 



Pip + 2) 



p = 4,6,8 



(6.3) 



(6.4) 



corresponding to the even unitary c < 3/2 minimal models [21|, [22|. Precisely this set of 
TBA equations had been conjectured earlier in [24jl, and then further generalized in [25 |. 

Finally, let us consider the model of primary interest here, namely, boundary SShG. For 
simplicity, we fix 0o = in the boundary Lagrangian ( |2.22| ), which corresponds to $ = 0.0 
Due to the roaming limit ( |6.3|) , we should rescale the remaining two parameters rj, (p so that 
the boundary entropy can be a function of well-defined (finite) boundary parameters. For 
this purpose we set r\ — > and ip — > oo while keeping OqI] and 6q — 2ip finite. Let us introduce 
new boundary parameters fi and f 2 defined by (see Eq. (|2.44|) ) 



29 oV 



7T 



h, l + -e^ 2 ^cosh/ 2 . 



We can reexpress the boundary entropy ( |5.23|) in terms of these parameters as 



SB 



+ S 



(2) 
B ) 



where 



— J d9^(9;f t ) L t (9), 



1,2, 



(6.5) 



(6.6) 



14 Consider the boundary SSG model first. When (j>o = the total Lagrangian respects C symmetry due 
to the Z2 symmetry (f> —* —<f>. Therefore, the boundary S matrix should respect C symmetry, namely the 
soliton and antisoliton should scatter equally on the boundary. Since the topological sector of the SSG S 
matrix is encoded in the SG part, the boundary parameter $ should vanish as it does in the SG model ||. 
This holds also for the boundary SShG S matrix because the two models are related by the fusion procedure. 
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with 

.„ „, 4 cosh 6 cosh f 

= cosh 20 + cosh 2/ ' (6 ' 7 > 

To compute the roaming boundary entropy, we fix a value of r where c e //(r) lies on a 
plateau (|6.4j) . Then, as we change the boundary roaming parameters f\ and f 2 , we check if 
the boundary entropy interpolates between the values [^7j 



s B {r,s) = In 



sm( — ) \ / sin l , 



= s B (r,l) + s B (l,s). (6.8) 

corresponding to conformal boundary states (r , s) (which, in turn, correspond to primary 
fields $( r , s )). 

Indeed, we can see clearly from Figure ^| |^] that Sg interpolates between boundary 



entropies of the conformal boundary states 

(l,a)^{ {a ~ 2 ' 1) , a odd . (6.9) 

Similarly s^' generates the new flow (see Figure H) 

f (a- 2,1) . . 

(l,a)^ J , a even . (6.10) 

While these flows are generated by changing one parameter while fixing the other, we can 

generate more general flows by changing fi and /a simultaneously. In view of the additivity 



15 Note that this expression satisfies Ss(l,l) = 0. The correct expression for the conformal boundary 
entropies has an additional "constant" term (i.e., independent of both r and s); we neglect this term here, 
since we are mostly interested in differences sb{t , s) — sb{t' , s'), for which the constant term cancels. 

16 We recall that for each bulk primary field < i ) (r,s)j there corresponds a conformal boundary state 
\h( r ,s)) (which, for brevity, we denote here by (r , s) ) such that the partition function Zn s ) for the CFT 
on a cylinder with conformal boundary states (1,1) and (r , s) is given by 

i.e., the character of $( r>s ). In particular, 

Z(i,i)(i,i) = X(x,x)(?) 

is the character of the unit operator. 

17 For p > 4, we cannot associate any conformal boundary state to the final plateau (i.e., for asymptotically 
large values of the boundary parameter / 2 ), since there is no state (0,1). 
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Figure 5: Boundary roaming trajectories: vs. f 



property ( |6.8|) , these two sets of flows can be combined to generate additional flows for the 
total boundary entropy sb 

(s,r) 

(s_2 ' r) , r-s= odd . (6.11) 
(s,r + 2) 

I (s-2,r + 2) 

Note that r — s= even/odd corresponds to the Neveu-Schwarz/Ramond sectors, respectively. 



(r,s) 



7 Discussion 

We have presented the exact solution of the boundary SShG model - an integrable QFT 
whose bulk and boundary S matrices are not diagonal. In particular, we have derived an 
exact inversion identity ( |4.5|) - ( |4.7|) , as well as the TBA equations and boundary entropy 
( [5.23Q . Moreover, we have uncovered a rich pattern of boundary roaming trajectories, which 
remain to be understood in detail. 

Although the boundary SShG model has a special feature which allows it to be solved by 
an inversion identity (namely, the bulk S matrix satisfies by free-Fermion condition (|2.11|) ), 
it is by no means the only such model. Indeed, there are infinite families of integrable QFTs 
with N = 1 or iV = 2 supersymmetry [[52]-|55| that have this property. These models have 



bulk and boundary S matrices which are similar to those of SShG, and therefore, we expect 
similar inversion identities to hold. We hope to report on these models in the near future 
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Figure 6: Boundary roaming trajectories: vs. f2- 



Finally, we recall [^6[ that one can readily obtain the Hamiltonian of an integrable open 
quantum spin chain with N spins from any homogeneous open-chain transfer matrix t(#|0) 
( |4.1| ). Indeed, the Hamiltonian 7i is given by 





H oc -t(0|O) 



(7.1) 



which commutes with t(6*|0). For the R matrices which we have considered here ( |2.9|) , (|2.30| ), 
the corresponding Hamiltonian is that of a certain anisotropic XY chain with both bulk and 
boundary magnetic fields. By determining the eigenvalues ( [4.17 ) of the transfer matrix, we 
have evidently also solved the corresponding open quantum spin chain. It would be inter- 
esting to exploit this solution to determine properties of this model in the thermodynamic 
limit. 
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A Relation of Yang matrix to Sklyanin transfer matrix 



In Section 3.2, we stated that the Yang matrix (|3.8|) is related to the Sklyanin open-chain 
transfer matrix (|3.9|) in the following way ( |3 . 1 1|) : 

YQ = T(6i\ei,... ,e N ), i = i,...,N. (a.i) 

We present here a proof for the case N = 2. Evaluating the transfer matrix at — 9±, we 
have 

r(0i|0i,0 2 ) = tr o {So(-0i+m^ + )* o 5o 2 (0i-0 2 )5oi(O)So(0i;e-)5oi(20i)5 O2 (0i + 2 )} 
= tro^oa^i - 2 )7>oiS o (0i ;e-)(^oi^Qi)^oi(20i)(7 3 oi^oi) 
x 5 O2 (0 1 + 2 )(PoiPoi)S o (-0 1 +z7r;e + ) to } = --- (A.2) 

In passing to the second line, we have used the cyclic property of the trace, as well as 
5(0) = V and V 2 = I, where V is the permutation matrix ( [2.14 ). 

• • • = S l (9 1 ; £_) tr o {5 O2 (0i - 2 )5 O i(20 1 )S' 12 (0 1 + 2 )P O iS o (-0i + «r ; £+)'«} = ■ ■ ■ (A.3) 
Here we have used P i A Poi — Ai, and the P symmetry of the R matrix fl2.13|) . 

• • • = S x (0i ; £_) tr o {S 12 (0! + 02)501(200502(01 - 2 )P O iS o (-0i + in ; £+)' } = • • • (A.4) 
Here we have used the Yang-Baxter equation ( |2.12j) . 

■ •• = Si(0i ;e-)5 12 (0i + 2 )tr o {5 O i (20i)(P O i^oi)5 O2 (0i - 2 )P O iS o (-0i + in ; £+)' } 

= Si(0i ;e-)5i 2 (0! + 2 )tr o {5 O i(20i)PoiSo(-0i + ivr ; e+)*°}5 12 (0 1 - 2 ) 

= Si(0i;e-)5i 2 (0i + 2 )Si(0i;e+)5i 2 (0i-0 2 ). (A.5) 

In passing to the last line, we have used the boundary cross-unitarity relation ( 2.34j) with 
= Y — 0i, and the crossing relation Soi(in — 0)* 1 = 5oi(0). Comparing the last line to the 
expression ( |3.8| ) for the Yang matrix, we conclude that 

r(0 1 |0 1 ,0 2 )=y (1) . (A.6) 

For higher values of N, the proof is similar. 



B Derivation of inversion identity 

In Section 4.1, we give the important inversion identity ( f4.5|) - ( |4.7| ). Here we explain in 
more detail how we derived it. As already mentioned in text, the main idea is to formulate 
the fusion formula, following Ref. p2|| , to which we shall refer as I. Q 
18 In order to facilitate comparison with |32j, we use here similar notations. 
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Although the "dressed" bulk S matrix S(9) (|2~20D is regular at 9 = 0, the "bare" bulk S 
matrix R(9) (|2.9| ) has a pole there. In order to avoid complications from this spurious pole, 
in this Appendix we rescale R(9) by the factor sinh#; i.e., we take R(9) to be given still by 
T9|), but now with matrix elements 



a ± {6) 
c{9) 



± sinh 9 — 2i sin Bit . 

Q 

—2i sin Bit cosh - , 



b{9) = sinh#, 
d(9) = —2 sin Bn sinh 



9 



(b.i; 



Keeping in mind the symmetries ( |2.13| ) of the R matrix, the unitarity relation (I 2.3) is 



Ru{9) R 12 {-9)=({9)I, 
and the crossing relation (I 2.4) is 

Rn(9) 

withQ 



e 



e 



((9) = -4 coslr -(sinlr - + sin 2 Bit 



V 1 R 12 (-9-p) t2 V 1 



i ■ 



P 



Z7l 



V 




(B.2) 



(B.3) 



(B.4) 



The matrix Ru(9) at 9 



—p is proportional to the one-dimensional projector P l2 
( 1 -1 \ 



12 



V 



ooo 


-10 






1 / 



12) 



p 



12 



(B.5) 



As explained in I, from the corresponding degeneration of the (boundary) Yang-Baxter equa- 
tion, one can derive identities which allow one to prove that fused (boundary) S matrices 
satisfy generalized (boundary) Yang-Baxter equations. 

The fused R matrix is given by (I 2.13) 

R<n>3(0) = P+ R l3 (9) R 23 (9 + p) P+ , (B.6) 

where P^ 2 = I — P{ 2 - An important observation (which one can verify by direct calcula- 
tion) is that the fused R matrix can be brought to upper triangular form by a similarity 
transformation 



X12 R<i2>s(9) Xu = upper triangular 



(B.7) 

19 Alternatively, choosing p — —iir, one has V — I. 

20 This observation is similar to, but not the same as, the one made by Felderhof [Q. Indeed, in our 
language, he shows that Ri3(0)R23(9 + p) (i.e., the expression for the fused transfer matrix without the 
projectors P-jt,) can be brought to triangular form by a (somewhat more complicated) ^-independent similarity 
transformation. Although for the case of periodic boundary conditions both approaches lead to the inversion 
identity, this appears to be no longer true for the case of boundaries. 
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where the 4x4 matrix X is independent of 6, and is given by 



X 



I* 






-sin^ cos^ 
cos ^ sin ^ 



B~ 



73 ^ 



\73 



X =1. 



:b.s) 



-73/ 



It follows that the fused monodromy matrices (I 4.7), (I 5.4), (I 5.5) 

T<:i2>(9) = R<12>n{9) • ■ ■ R<i2>i{9) ; 
^<12>(# + P) = -R<i2>i(#) • • • R<12>n{@) , 



(B.9) 



also become triangular by the same transformation. 

Denoting (as in I) our "bare" boundary S matrices R(9 ; </?_), R(— + ; f + ) by K~(6), 
K + (9), respectively, the corresponding fused matrices are given by (I 3.5), (I 3.9) 



^<i 2 >(0) = P? 2 K- (6) R 12 (26 + p) K 2 (6 + p) P+ , 
K+ 12> (9) = {P+ K+iOf 1 Rui-20 - 3p) K+(9 + p) t2 P+} tl2 



;b.io) 



since M = F*V = I. 



Remarkably, the fused K matrices are also brought to upper triangular form by the same 
similarity transformation 

X12 K^ 12> {6) X\2 = upper triangular . (B.ll) 

It follows that the fused transfer matrix t(9), which is given by (I 4.5), (I 4.6) 

1(6) = tr 12 K+ 12> (9) T <12> (9) K~ 12> (9) T <12> (9 + p) , (B.12) 

is proportional to the identity matrix, 

t(0)ocl, (B.13) 

where the proportionality factor is determined from the diagonal elements of the various 
triangular matrices. 

The fusion formula is given by (I 4.17), (I 5.1) 

1 



t(9) i{9 + p) 



C(2# + 2p) 



For simplicity, we consider here the homogeneous case 



t(6) + A{K + (6)}A{K-(6)}5{T(6)}5{f(6)}] , (B.14) 

, AO- 



ID, i=l 
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where the transfer matrix t(9) is given by (|4.1| ) (see also (I 4.1), (I 4.2)), and the quantum 
determinants |50| are given by (I 4.15), (I 5.3), (I 5.7) 



5{T{9)} = 5{f(6)}=a6 + pf, 

A {K-{9)} = tr 12 {Pf 2 K^{9) R 12 (29 + p) K^{9 + p) V x V 2 } , 

A {K+{9)} = tr 12 {Pf 2 V 1 V 2 K+{9 + p) R 12 {-29 - 3p) K+{9)} . (B.15) 

Reverting to the original normalization of the R matrix by rescaling each of the transfer 
matrices t(9) in ( [B. 14| ) by (sinh^) -2 ^, introducing the inhomogeneities 9i in the obvious 



way, and factoring the result into a product of two factors, we arrive at the results ( |4.5|) - 
(0). 
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